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Abstract: The new Rain-Sea Interaction Facility (RSIF) has been established at GSFC/WFF 
and the first findings are presented in this report. The unique feature of this laboratory is the 
ability to systematically study microwave scattering from a water surface roughened by 
artificial rain, for which the droplets are at terminal velocity. The fundamental instruments 
and systems (e.g. the rain simulator, scatterometers and surface elevation probes) were 
installed and evaluated during these first experiments - so the majority of the data were 
obtained with the rain simulator at 1 m above the water tank. From these initial experiments, 
three new models have been proposed: the square-root function for NCS vs R, the log- 
gaussian model for ring-wave elevation frequency spectrum, and the Erlang probability 
density distribution for backscattered power. Rain rate is the main input for these models, 
although the coefficients may be dependent upon other factors (drop-size distribution, fall 
velocity, radar configuration, etc.). The facility is functional and we foresee collaborative 
studies with investigators who are engaged in measuring and modeling rain-sea interaction 
processes. 


Avant-propos: Une experimentation de pluie simulee a ete faite pour etendre a forts taux de 
pluie (R mm hr' 1 ) le modele a faibles taux propose par Bliven pour la section radar efficace 
normalisee (NCS). La dependance de la NCS est lineaire en fonction de faibles taux de pluie. 
II a ete demontre que pour ces conditions d'experimentation la theorie de retrodiffusion de 
Bragg pour les vagues annulaires est applicable. Ce papier rectifiera le domaine de validite du 
modele lineaire et proposera un modele plus general applicable aussi bien pour faible pluie 
que pour forte pluie. 

En vue d'accomplir cette tache, on proposera un modele original pour la distribution 
de densite de probability des donnees radar. Egalement une etude de la topographie de surface 
est faite a l'aide d'une sonde d'elevation par capacitance, qui fournira les donnees necessaires 
pour une interessante moderation de spectre d'elevation de surface en frequence dans 
I'intervalle d'etude des vagues de type capillaire-gravites. 

Le simulateur de pluie qui a servi a mener cette experience sera decrit pour le lecteur 
dans le but de lui permettre de decouvrir d'une part I'originalite de cette realisation et d'autre 
part la capacite a simuler au mieux une precipitation naturelle. 
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1. Introduction 

Rain and wind measurements in oceanic regions will make significant contributions to 
weather monitoring and climate studies. Spaceborne instruments such as altimeters and 
scatterometers measure power levels of microwave echoes from the sea surface and since 
echo strength correlates with sea-surface roughness, particularly small-scale roughness, any 
processes that affect sea surface roughness can modify measurements. Nearly ubiquitous 
winds roughen the sea surface, so most data inversion algorithms deal with wind exclusively. 
Rain agitates the sea surface. Atmospheric attenuation and scattering by rain are included in 
two models that simulate satellite-based systems (at nadir by Meneghini and Atlas (1986) and 
at scatterometer angles by Sobieski et al. (1991)) - but the effects of rain on microwave 
scattering from the sea surface are yet to be included in space applications. 

When a raindrop hits a water surface, it can generate a cavity with a crown, which 
collapses to form a vertical stalk of water, which subsides to spawn rings of gravity-capillary 
waves that propagate outward (Worthington (1882), Le Mehaute et al. (1987) and Le Mehaute 
(1988)). At grazing angles, Wetzel (1990) found from analysis of laboratory data that stalks 
are the dominant feature contributing to backscattered power. At incidence angles used from 
space, a dearth of data hampers validation of numerical models. Two recent studies by Bliven 
and Giovanangeli (1993) and Bliven et al. (1993) use Ku- and Ka-band scatterometers (30 
degrees incidence angles and vertical polarization) to study scattering from a rain-roughened 
water surface for simulated light rains (R < 30 mm hr' 1 ) in a wind-wave tank. Those data 
show (a) that the dependence of Normalized Radar Cross-Section (NCS) on rainfall rate R is 
well modeled by a linear function and that (b) there is evidence that the dominant scattering 
mechanism is Bragg scattering from ring-waves. Much more can be learned about microwave 
scattering from rain-roughened surfaces by conducting experiments that span a broader range 
of physical conditions. Higher rain rates obviously need to be studied and other rain-related 
topics include the effects of dropsize and fall velocity. On the other hand, microwave topics 
include scattering physics related to radar viewing angle, polarization, and microwave 
frequency. 

A research facility was needed to allow these issues to be systematically investigated 
using both standard and specialized instrumentation in a setting that permits water drops to 
reach terminal velocity. Such a site was available at the Wallops Flight Facility and as part of 
this endeavor, the Rain-Sea Interaction Facility was established. Research being conducted at 
the RSIF is part of an ongoing international research program that includes numerical 

modeling of microwave scattering from the sea by scatterometers and altimeters (Sobieski et 
al. 1990) and simulations of scattering from rain- and wind-roughened seas in a laboratory 
wind-wave tank (Giovanangeli et al. 1991). The new RSIF is functional and open for 
cooperative investigations. 

In this paper, we present a guide to papers pertaining to rain-sea interaction processes, 
a description of the new Rain-Sea Interaction Facility, and the first data sets from the RSIF 

are examined. The main objective was to obtain data over an extended range of rainfall rates 
to see if the linear relationship between NCS and R is robust. Two additional topics are 
addressed: characterization of ring-wave frequency spectra as a function of R and the 
functional form of backscattered signal probability-density-distribution as a function of R. 
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2. Background 

Case studies show that rain effects data from altimeters, scatterometers and SAR's 
(Guymer et al. 1981, Fu and Holt 1982, Hawkins and Black 1983, and Black et al. 1985, 
Srokosz and Guymer 1988) and that atmospheric corrections do not necessarily remove all 
biases (Black et al. 1985). So air-sea interaction processes need to be investigated. Laboratory 
data sets (Moore et al. 1979 and Bliven et al. 1989) indicate that (a) Ku- and Ka-band 
scatterometer backscattered power can be enhanced by water drops striking the air-water 
interface, and (b) that water surfaces agitated by solely light rains can yield about the same 
backscattered power as water surfaces roughened by solely light winds. Moreover, for 
combined light winds and light rains, Bliven and Giovanangeli (1992) established that 
backscattered power can be well modeled as the sum of rain and wind contributions. At 
higher winds and rain rates, various physical processes may be important. Attenuation of wind- 
waves by rain-generated turbulence in the water has been measured by Tsimplis (1992), who 
found that estimates of eddy viscosity coefficients by Manton (1973) are too large by two 
orders of magnitude, so it is doubtful that wave-turbulence interaction is the dominant 
mechanism for wind-wave attenuation by rain. On the other hand, a model by Le Mehaute 
and Khangaonkar (1990) indicates that raindrop impacts on the water surface can be 
inhomogeneous and consequently, rain may contribute to long-wave growth. Gravity-wave 
interaction with capillary-waves has been modeled by Kharif et al. (1990), who report 
asymmetries as a function of the angle between the propagation directions of the two wave- 
trains. For application to in situ conditions, this model needs to be extended to include the 
effects of randomly distributed drop impacts. Consequently, surface agitation by rain needs to 
be included in empirical data inversion algorithms and numerical models — but data are 
certainly necessary to guide development and assessment of theoretical models. 

Since the pioneering investigations by Wright (1966 and 1968) and Bass et al. (1968), 
estimates of radar backscattered power for viewing angles in the range of 30 to 60 from nadir 
have been typically computed by representing the sea surface as a two-scale rough surface. 
The strengths and limitations of this formulation for wind-wave studies were reviewed by 
Plant (1990), who discusses that composite models are used to represent wind-roughened seas 
as patches of short waves (Bragg scattering model) that are locally tilted by large wave 
components. There is considerable evidence that scatterometer measurements are highly 
correlated to wind friction velocity (Jones and Schroeder 1977, Li et al. 1989, Giovanangeli et 
al. 1991 and Bliven et al. 1992). There also is considerable evidence that short gravity and 
capillary waves are the dominant spectral region that supports the bulk of the momentum 
fluxes from the wind to the water surface (see Plant (1982) for a review and Geernaert et al. 
(1986) for a recent assessment). Both aircraft data obtained by Ross and Jones (1978) and 
SEASAT data analyzed by Glazman (1988) show that the effect of long waves on 
scatterometer average backscattered power is small. With a recognition of these factors, the 
entire wavenumber spectrum is not needed for first-order modeling of average quantities. 

Indeed, Keller et al. (1992) employed laboratory results to assess assumptions concerning 
physical processes used in scatterometer wind-retrieval algorithms developed for field 
applications. Data from long fetches in wide wind-wave tanks can assist in development of 
application algorithms for combined wind and rain conditions. 
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The Bragg scattering theory should be useful for modeling scattering from ring-waves. 
For example, scattering features generated by a single drop impacting a calm water surface 
were investigated by Wetzel (1990), who concludes that due to its small amplitude, "a 

perturbation approximation may legitimately be applied at microwave frequencies," and who 
uses Wright's (1966) model to effectively replicate backscattered power time series from a 
ring-wave. Although Bragg scattering of microwaves from short wind-waves is widely used as 
the dominant mechanism in scatterometer wind algorithms and numerical models, the role of 
Bragg scattering from ring-waves has been quantified only for a restricted set of conditions by 
Bliven et o/.(1993). A complete survey of rain effects on sea-surface processes is beyond the 
scope of this report, however, an extensive list of references is provided in section 8. Rain- 
generated cavities, stalks, and ring-waves need to be measured simultaneously with 
microwave measurements in order to evaluate the role of these features. 

3. RSIF description and experimental conditions 

The Rain-Sea Interaction Facility at WFF was created during the summer of 1993 in 
order to make measurements of microwave scattering from water surfaces agitated by 
artificial rain. The purpose of this facility is to provide a controlled environment to study rain 
effects on microwave scattering from water surfaces and to use this data to develop and 
improve remote sensing algorithms. The aircraft hangar known as building N159 (Figure 3.1 ) 
contains a tower that was originally used for drying parachutes. Its dimensions are 4 x 4 x 17 
m and there is a water drain located in the center of the floor. A ladder in one comer leads to 
the high cat-walk platform at 14 m. This area was ideally suited for conversion to the Rain- 
Sea Interaction Facility. 

We installed equipment at ground level for these studies as shown in Figure 3.2. The 
water tank (2 x 1 x 1 m) is centered above the drain and for operations it is filled to a depth 
of 80 cm. Prior to each experiment, the tank is cleaned and filled with fresh water. A 6-cm 
hose is connected to the drain at the bottom of the tank and it is extended upward for use as a 
water overflow outlet, so that the water level is maintained at a constant height during 
experiments. The 13.5-GHz and 36-GHz scatterometer transmitter and receiver horns are 
pointing towards the center of the water tank. The storage shelves house the radar systems 
electronics, the capacitance probe electronics and a PC-based data acquisition system. The 
scatterometer systems are calibrated relative to the backscattered power from a 15-cm metal 
sphere that is located at the operational range to the water surface (1 m). The water tank is 
empty during this calibration procedure. A capacitance wire-probe is used to measure water 
surface time-series in the tank at a point adjacent to the rain-simulator footprint. Data quality 
is assured by analyzing characteristics of each data set during each experiment. 
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The high-capacity rain simulator ( Figure 3.3) was built so that a broad range of rain 
rates could be simulated, especially extreme rainfalls. It consists of a stainless steel box (80 x 
80 x 6 cm) that has a water inlet and an air-vent that is open during filling — but closed 
during operation so that the system is sealed. A pump with volume control regulates the flow 
rate of filtered water to the rain simulator. The bottom of the box ( Figure 3.4) contains 1100 
holes spaced with 2.5 cm between center points. A plastic nipple is inserted into each hole 
and to these, various sized hypodermic needles can be attached from the exterior. Water drops 
with diameters ranging from 1.2 to 4.2 mm diameter are produced by using various sized 
needles. The rain simulator is centered above the water tank and it can be lifted and attached 
to the cat-walk by means of a 500- lb-capacity winch that is attached to the cat-walk. 

For this investigation, analog signals from the radars and capacitance probe were 
digitized at 256 Hz for 128 s to obtain 32k of data for each experimental condition. 

4. NCS versus R: Square-root model 

Radar time-series of backscattered power measurements from the water surface 
roughened by 2.8-mm diameter drops falling from 1 m are shown in Figure 4.1. The relative 
power RP values are normalized by the return power level from the calibration sphere. For 
similar conditions, the backscattered power levels are usually lower for the 36-GHz system 
than for the 13.5-GHz system. 

Average backscattered power levels NCS have been computed from the times-series 
data and a typical summary figure for the 13.5-and 36-GHz scatterometers is shown as a 
function of R for rain rates up to 200 mm hr" 1 in Figure 4.2. These results were obtained 
using 2.8-mm-diameter drops falling from 1 m. To gain experience and develop techniques, 
we conducted numerous experiments with the rain simulator at 1 m from the calm water 

surface. One data set was obtained with the simulator at 13 m, but that series was terminated 
in order to make modifications to the simulator when a seal suddenly burst. Analysis of the 
data from 1 m show trends that are similar to those from the data from 13 m. Certainly 
investigation of drop-size and fall velocity effects are appropriate and will be conducted in the 
future — but much can be gained from analyzing the 1-m data and we present those results 
which are a guide to future studies. For the square-root model, a small DC-offset has been 
removed to force the model to zero when there is no rain. The data are well described by a 
square-root function, as shown in Figure 4.2 for 2.8-mm-diameter drops falling from 1 m. 

At light rainfall rates, the square-root relationship is similar to a linear model — but at 
higher rainfall rates, the relative increase of NCS with R decreases so a linear model is not 
appropriate. For the entire range of simulated rains, a square-root model represents the Ku- 
and Ka-band data obtained with the rain simulator (2.8-and 4.2-mm-diameter drops) at 1 m 
from the water surface very well. Preliminary results with the rain simulator (2.8-mm drops) 
at 13 m exhibit a similar trend. It is interesting to observe that at high rainfall rates, NCS 
continues to increase as R increases — i.e., there is no evidence of saturation. Thus all of these 
data from an extended range of rainfalls indicate that a square-root relationship between NCS 
and R is fitting. 


4 



5. Ring-wave spectra: Log-Gaussian model 

To obtain a quantitative measurement of the rain-roughened water surface, we used a 
capacitance probe to measure surface elevation time-series at a point adjacent to the rain 
simulator footprint. A typical time series is shown in Figure 5.1 which indicates that although 
the ring-waves are small (typically less than 1 mm in height), the measurement is not noisy. 

A similar arrangement was used by Bliven et al. (1993b), however, that rain simulator 
(also about 75 x 75 cm 2 ) contained only 72 nozzles and the capacitance probe was located 
exceedingly close to the impact point of drops from one of the nozzles. As a result, the 
surface elevation frequency-spectra contained features related to the drop impact-rate (number 
of drops per second striking the surface) as well as the more random ring-waves from the 
other 71 nozzles. Unfortunately even light rains required a drop production rate for each 
nozzle that was higher than 2 Hz so the two spectral features are not easily separable — thus 
interpretation of those data is hampered. One of the design specifications for the new rain 
simulator at RSIF was that rainfall rates up to 200 mm hr' 1 be feasible with 2.8-mm-diameter 
drops at an average production rate less than 2 Hz from the nozzles. This specification 
facilities the unambiguous identification of ring-wave components in surface-elevation 
frequency spectra. 

Typical ring-wave frequency spectra are displayed in Figure 5.2, from which it is 
apparent that the peak frequency is about 6.5 Hz, the bandwidth is quite narrow and the decay 
is asymmetric — faster decay on the low-frequency side than on the high-frequency tail. 

We model these spectral features by using a logarithmic gaussian relationship given 
by: 
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where: 

S,,, : spectral peak magnitude. 

f p : peak frequency. 

Af : frequency bandwidth. 
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These three parameters can be obtained from the data by a least-square algorithm, examples 
of the close agreement between the model and the spectral data are shown in Figure 5.2. This 
analysis has been done for all the rainfall-rates and the results are shown in Figure 5.3, which 
shows that the spectral peak frequency is a weak function of rainfall-rate and the spectral 
bandwidth is approximately constant. The 3-dB bandwidth and the bandwidth parameter Af 
are related by 


A 4/5 


2 f p sinh 


' Af 

ln(2) ' 

<P \ 

n 


( 2 ) 


In spite of that complicated relationship between Af and Af MB , they are approximately equal 
(Af is ~5% greater than the 3-dB bandwidth). Thus another feature of this model is that Af.S fp 
contains more than 90% of the total energy, so a close approximation of the elevation 
variance can be estimated from two of the characteristic parameters, i.e. 


£ - 3, a / 


(3) 


where the subscript m means relative to the model. 
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Wavenumber elevation spectra can be obtained from frequency spectra by using an 
appropriate dispersion relationship which depends upon the fluid properties and wave 
characteristics. For gravity-capillary waves, a commonly used dispersion relationship is 


a 


2 



k 3 


where: 


k : wavenumber spectrum 

o=2nf : radian frequency 

g : gravitational constant 

p : water density 

x : water surface tension 


( cm ‘ ) 

( rd s' ) 

(980 cm s' 2 ) 

( 1 g cm 3 ) 

(~74 dynes cm' 1 ) 


(4) 


The change of variables between frequency space and wavenumber space is accomplished by 
using a dispersion relationship and the following formula 
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where the derivative of radian frequency with respect to wavenumber is the phase speed V. 

Using equation (4), the log-gaussian frequency spectral model can be expressed in 
wavenumber space as 
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Although this analytical wavenumber spectral model is derived directly from a simple 
frequency spectral model, it is cumbersome to use, particularly because the k dependence 
appears in the amplitude factor V(k). 
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A simpler spectral model can be written as 
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with the coefficients given by 
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( 10 ) 


where V p is the phase speed at the peak frequency. 

This simple model provides a good approximation of equation (6) for k ranging from 
0.04 rd cm 1 to 4.3 rd cm' 1 , i.e. it has a similar shape and represents 98% of the total energy. 

It is interesting to note that for these experimental data, f p is -6.4 Hz so V p is close to the 
minimum phase speed V„ (17.83 cm s.j and k p -1.43 rd cm' 1 , which is consistent with a 
model developed by Le Mehaute (1988) for individual drops and extended by Sobieski (1991) 
for natural conditions. 

A typical surface elevation probability density distribution is shown in Figure 5.4, and 
analysis of all the surface elevation pd indicates that the distributions are close to Gaussian. 

The elevation psd for 11.6 Hz was used to estimate the psd for the Bragg resonant 
wavelength for the 13.5-GHz scatterometer. These psd values are shown in Figure 55, and a 
square-root dependence on R represents the data well. We note that both the elevation 
variance and elevation psd at 11.6 Hz have trends which are consistent with the NCS square- 
root vs R model. 
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6. RP probability density distributions: Erlang model 


The scatterometer RP time-series (r) are positive signals ( Figure 6.1) and the 
probability density distributions look like the familiar Rayleigh distribution ( Figure 6.2). 
Relative power RP values are normalized by the return power level from the calibration 
sphere. Each probability density distribution increases from zero to a maximum value and 
decrease toward zero as r increases. Unfortunately a semi-logarithm plot shows that the decay 
is linear in log scale and not parabolic as it is for Rayleigh function. So the decay is 
exponential with respect to r higher than the peak value in regular scale. The Rayleigh 
distribution is totally determined by a single parameter — either the mean value or the standard 

deviation — so it is impossible to independently adjust the magnitude of the peak and the 
decay rate. 

The Erlang distribution is a two-parameter function so it can account for a linear 
decay in log scale and variable magnitude of the peak. The form of the Erlang distribution is 




p -Xr 

r(a+l) 


(ID 


where: 

r : normalized radar data with respect to the sphere, 
a: order of Erlang function. 
k: exponential coefficient. 

T : Euler gamma function or generalized factorial function. 

The first and second moments of the distribution, i.e. the mean <r> and the standard 
deviation a „ are functions of the two Erlang's parameters. 


mean value: 


standard deviation: 


Likewise, a simple formula relates alpha 
inversion of the above equations. 
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and lambda to mean and standard deviation 


by 


Alpha parameter : 
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(14) 


Lambda parameter: 



(15) 
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The mean values and standard deviations of RP are well modeled by a square-root 
function over a broad range of rain-rates (10 to 210 mm hr *) as shown in Figure 6.2. A 
small DC offset has been removed from the square-root model to force the model to zero 
when there is no rain. From the RP data time-series, mean and standard deviation values were 
computed and used to model the measured PDs by the Erlang distribution. The close 
agreement with the observations is apparent in Figure 6.3. 

For both radars, alpha is an increasing function. On the other hand, lambda is a 
decreasing function of rain-rate. This combination of parameters leads to the following 
observations. High-RP values ( 'big' features) have higher probability to occur as rainfall-rate 
increases, therefore rare events at light rain become common events at heavy rain. 

7. Summary 

The new Rain-Sea Interaction Facility has been established and the first finds are 
presented in this report. The unique feature of this laboratory is the ability to systematically 
study microwave scattering from a water surface roughened by artificial rain, for which the 
droplets are at terminal velocity. The laboratory instruments and systems such as the rain 
simulator were evaluated for functionality during these first experiments so the majority of the 
data were obtained with the rain simulator at 1 m above the water tank. The tests were 
successful, so future experiments will be able to use the rain simulator attached to the cat-walk 
at a height of 13 m. Preliminary data from this evaluation period indicate that the data trends 
from 1 m do not change dramatically when the rain simulator is at the full height. Data are 
needed, however, to develop/validate scaling laws for non-terminal velocity conditions, which 
are typical for laboratory settings. The log-gaussian spectral model and the Erlang probability 
density distribution model should be able to represent data from a broad range of conditions, 
nonetheless, systematic studies need to be conducted to obtain parameter coefficients from 
configurations closer to natural conditions. Thus the results presented here are intended to 
illustrate the capability of the new Rain-Sea Interaction Facility and further investigations are 
planned to provide improved understanding of physical processes. 

From these initial experiments, three new models have been proposed: the square-root 
function for NCS vs R, the log-gaussian model for ring-wave elevation frequency spectrum, 
and the Erlang probability density distribution for backscattered power. Rain rate is the main 
input for these models, although the coefficients may be dependent upon other factors (drop- 
size distribution and fall velocity). Thus further experiments are needed to refine the models. 
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146-70-11) contributed funding for this research. Jean P. Schermann of the Formation 
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Figure 1. Rain-Sea Interaction Facility (17 m white tower). 
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Figure 2. Interior of Rain-Sea Interaction Facility: Radars and Tank. 
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Figure 4. Rain-simulator detail - 1100 nozzles with 2.54 cm center point spacing. 
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